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' Abstract 

Elek and Lippner (2010) showed that the convergence of a sequence of bounded-degree 
graphs imphes the existence of a hmit for the proportion of vertices covered by a maximum 
matching. We provide a characterization of the hmiting parameter via a local recursion 
' defined directly on the limit of the graph sequence. Interestingly, the recursion may admit 

Qh . multiple solutions, implying non-trivial long-range dependencies between the covered ver- 

r~| ' tices. We overcome this lack of correlation decay by introducing a perturbative parameter 

(temperature), which we let progressively go to zero. This allows us to uniquely identify the 
correct solution. In the important case where the graph limit is a unimodular Galton- Watson 
tree, the recursion simplifies into a distributional equation that can be solved explicitly, lead- 
ing to a new asymptotic formula that considerably extends the well-known one by Karp and 
04 , Sipser for Erdos-Renyi random graphs. 

> 

(N 

! 1 Introduction 

p 

I A matching on a finite graph G = {V, E) is a subset of pairwise non-adjacent edges M <Z E. 

The \V\ - 2\M\ isolated vertices of {V, M) are said to be exposed by M. We let M{G) denote 
the set of all possible matchings on G. The matching number of G is defined as 



KG) = ^max |M|, (1) 

and those M which achieve this maximum - or equivalently, have the fewest exposed vertices - 
are called maximum matchings. The normalized matching number of G is simply i/(G)/|y|. 

Our results belong to the theory of convergent graph sequences. Convergence of bounded 
degree graph sequences was defined by Benjamin! and Schramm [7], Aldous and Steele 12], see 
also Aldous and Lyons |lj . The notion of local weak convergence has then inspired a lot of work 
[H]) [H]) [E], [IS], [25]. In [17|, it is shown that the convergence of a sequence of bounded-degree 
graphs guarantees the existence of a limit for their normalized matching numbers. We provide 
a characterization of the limiting parameter via a local recursion defined directly on the limit of 
the graph sequence. In the important case where the graph limit is a unimodular Galton- Watson 
tree, the recursion simplifies into a distributional equation that can be solved explicitly, leading 
to a new asymptotic formula. 
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A classical example in this context is the Erdos-Renyi random graph with average degree c 
on n vertices, denoted by G{n,c/n) : as n — )• oo, G{n,c/n) converges in the local weak sense 
to a Galton- Watson tree with degree distribution Poisson with parameter c. In this case, Karp 
and Sipser ^2] showed that 

iy{G{n,c/n)) ^ ^ _ tg + e"'^*'^ + ctce~^*'= 



n n— s-oo / 

where tc £ (0, 1) is the smallest root of t = e^^^ (we will see in the sequel that the convergence 
is almost sure) . The explicit formula ([2|) rests on the analysis of a heuristic algorithm now called 
Karp-Sipser algorithm. The latter is based on the following observation : if e G -E is a pendant 
edge (i.e. an edge incident to a vertex of degree one) in G = (y,E), then there is always a 
maximum matching that contains e, so all edges that are adjacent to e may be deleted without 
affecting y{G). The first stage of the algorithm consists in iterating this until no more pendant 
edge is present. This is the leaf-removal process. G is thus simplified into a sub-graph with 
only isolated vertices, matched pairs, and a so-called core with minimum degree at least 2. As 
long as that core is non-empty, one of its edges is selected uniformly at random, the adjacent 
edges are deleted, and the whole process starts again. When the algorithm stops, the remaining 
edges clearly form a matching on G, but its size may be far below fiG) due to the sub-optimal 
removals on the core. 

On G(n, c/n), the dynamics of the deletion process can be approximated in the n — )• oo limit 
by differential equations which can be explicitly solved. In particular, the asymptotic size of 
both the optimal part constructed in the first stage, and the sub-optimal part constructed on 
the core can be evaluated up to an o(n) correcting term (which has been later refined, see [3]). 
Moreover, the second part happens to be almost perfect, in the sense that only o(n) vertices are 
exposed in the core. This guarantees that the overall construction is asymptotically optimal, and 
the asymptotic formula for v{Gn) follows. More recently, the same technique has been applied 
to another class of random graphs with a fixed log-concave degree profile [lOj, resulting in the 
asymptotical existence of an almost perfect matching on these graphs : 

I Vn I n^oo I 

In both cases, the proof of optimality - and hence the asymptotic formula for z^(G„) - relies on 
the fact that the second stage exposes only o(n) vertices, which is bound to fail as soon as one 
considers more general graph ensembles where the core does not necessarily admits an almost- 
perfect matching. We give simple examples in the Appendix. By using a completely different 
approach - namely establishing and solving an appropriate recursive distributional equation (a 
usual ingredient of the objective method, see [3]) -, we manage to obtain a general formula that 
considerably extends the above results. 

The rest of our paper is organized as follows: we state our main results in Section [2j In 
Section [3l we extend the Boltzmann-Gibbs distribution over matchings on a finite graph to 
infinite graphs. This will allow us to derive our Theorem[T]in Section^ We deal with the specific 
cases of trees (and random graphs) in Section [5l We end the paper with an Appendix presenting 
simple examples of graphs for which the limiting local recursion admit multiple solutions. 
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2 Results 



Let us start with a brief recall on local weak convergence (see [312] for details). A rooted graph 
(G, o) is a graph G = {V, E) together with the specification of a particular vertex o , called 
the root. We let Q denote the set of all locally finite connected rooted graphs considered up to 
rooted isomorphism, i.e. (Gi, oi) = (G2, 02) if there exists a bijection 7: Vi ^ V2 that preserves 
roots (7(01) = 02) and adjacency {uv G Ei <;=^ 'y{u)j{v) E E2). In the space G, a sequence 
{{Gn,°n)'i n £ 'N} converges locally to (G, o) if for every radius k £ 'N, there is € N such that 

n > Uk [Gn,On]k = [G,o]k. 

Here, [G, o]^. denotes the finite rooted subgraph induced by the vertices lying at graph-distance 
at most k from o. It is not hard to construct a distance which metrizes this notion of convergence 
and turns G into a complete separable metric space. We can thus import the usual machinery of 
weak convergence of probability measures on Polish spaces (see e.g. [9j). We define 7^(G) as the 
set of probability measure on Q. There is a natural procedure for turning a finite deterministic 
graph G = {V, E) into a random element of Q : one simply chooses uniformly at random a vertex 
o G y to be the root, and then restrains G to the connected component of o. The resulting law 
is denoted by h({G) G ^{Q). If (G„)nGN is a sequence of finite graphs such that {U{Gn))n&] 
admits a weak limit p G 'P{Q), we call p the random weak limit of the sequence {Gn)n&n- Finally, 
for any (i > 0, we define Qd as the space of all rooted connected graphs with maximal degree no 
more than d. 

Rather than just graphs G = {V,E), it will be sometimes convenient to work with discrete 
networks G = {V, E, Ai), in which the additional specification of a mark map : -E — >• N allows 
to attach useful local information to edges, such as their absence/presence in a certain matching. 
We then simply require the isomorphisms in the above definition to preserve these marks. 

The first main implication of our work is that the local weak convergence of a sequence 
of graphs is enough to guarantee the convergence of their normalized matching numbers to a 
quantity that can be described directly on the random weak limit of the graph sequence. 

Theorem 1 Let Gn = (Vn,En),n gN, be a sequence of finite graphs admitting a random weak 
limit p. Then, 

I ^ n I n— >oo 

where 7 G [0, ^] is defined by a recursion defined directly on the random weak limit p. 

Since the work of Heilmann and Lieb [2T], it is known that the thermodynamic limit for 
monomer-dimer systems exists and basic properties of this limit are derived for lattices. In 
particular, [211 Lemma 8.7] shows the convergence of the normalized matching number when 
the underlying graph is a lattice. More recently, Elek and Lippner [T7j extended this result by 
using the framework of local weak convergence for bounded degree graphs. Here we remove 
the bounded degree assumption. More importantly, the approach in [T7] is non-constructive. 
In contrast, we provide a characterization of 7 in terms of a local recursion defined directly 
on the random weak limit p. We postpone the discussion on how to actually compute 7 from 
p to Subsection 13.11 Our approach starts as in j21) with the introduction of a natural family 
of probability distributions on the set of matchings parametrized by a single parameter z > 
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called the Boltzmann-Gibbs distribution. The analysis in [2T] concentrates on the properties 
of the partition function, also known as the matching polynomial, from which a result like (jH) 
can be deduced. Our analysis differs from this approach and concentrates on the analysis of the 
local marginals of the Boltzmann-Gibbs distribution, in a similar spirit as in the (non-rigorous) 
work of Zdeborova and Mezard [26j. As in [26J and l_2D], we start from an elementary formal 
recursion satisfied by the matching polynomials, and deduce an exact recursion for the local 
marginals of the Boltzmann-Gibbs distribution on any finite graph. A careful analysis of the 
contractivity properties of this recursion allows us to define the monomer-dimer model on infinite 
graphs (see Theorem [6]), and to define their "normalized matching number" (see Theorem llOp . 
We should stress that the analysis of the marginal probabilities is essential for our second main 
result, namely the explicit computation of the matching number when the local weak limit is a 
Galton- Watson tree. Although simple adaptations of the argument in [21] would yield a result 
like Theorem [H the limit would be given in an implicit way which would not be sufficient to get 
our second main result. 

As many other classical graph sequences, Erdos-Renyi graphs and random graphs with a 
prescribed degree profile admit almost surely a particularly simple random weak limit, namely 
a unimodular Galton- Watson (UGW) tree (see Example 1.1 in [T|). This random rooted tree 
is parametrized by a probability distribution vr G 'P(N) with finite mean, called its degree dis- 
tribution. It is obtained by a Galton- Watson branching process where the root has offspring 
distribution tt and all other genitors have offspring distribution t: S 7^(N) defined by 

Vn £ N, TTn = ^ ■ . (5) 

Thanks to the markovian nature of the branching process, the recursion defining 7 simplifies 
into a recursive distributional equation., which has been explicitly solved by the authors in a 
different context |12j. 

Theorem 2 With the notation of TheoremUl if the random weak limit p is a UGW tree with 
degree distribution tt, we have the explicit formula 

1 - maxtg[o,i] F{t) 
^ = 2 ' 

where 

F{t) = t^'il -t) + <p{l-t) + <p(^l- ^^^^) - 1, 
and (j){t) = vr^t" is the moment generating function of the degree distribution n. 



Differentiating the above expression, we see that any t achieving the maximum must satisfy 

For Erdos-Renyi random graphs with connectivity c, the degree of the limiting UGW tree is 
Poisson with parameter c (i.e. (j){t) = exp(ct — c)), so that ([6]) becomes t = e~'^^ " . We thus 
recover precisely Karp and Sipser's formula ([2]). Similarly, for random graphs with a prescribed 



4 



degree sequence, the log-concave assumption made by Bohmann and Frieze guarantees that the 
above maximum is achieved at t = with F{0) = 0, hence ([3]) fohows automatically. 

A classical area of combinatorial optimization is formed by bipartite matching [23] . We end 
this section, with a specialization of our results to bipartite graphs G = {V = V"" UV^,E). The 
natural limit for a sequence of bipartite graphs is the following hierarchal Galton- Watson tree 
parameterized by two distributions on N with finite first moment, vr" and vr'' and a parameter 
A G [0, 1]. We denote tt" and tt^ the corresponding distributions given by the transformation ([5]). 
We also denote (p"" and c/)^ the generating functions of vr'^ and vr^. The hierarchal Galton- Watson 
tree is then defined as follows: with probability A, the root has offspring distribution tt"", all odd 
generation genitors have offspring distribution vr^ and all even generation genitors have offspring 
distribution vf^; similarly with probability 1 — A, the root has offspring distribution vr^, all odd 
generation genitors have offspring distribution vr" and all even generation genitors have offspring 
distribution vr^. In the first (resp. second) case, we say that the root and all even generations 
are of type a (resp. b) and all the odd generations are of type b (resp. a). To get a unimodular 
hierarchal Galton-Watson (UHGW) tree with degree distributions vr" and vr'', we need to have: 
\(t)^'{l) = (1 - A)(/>'"(1), so that 

,^ - 4>"{l) 



Theorem 3 With the notation of TheoremUl assume that the random weak limit p is a UHGW 
tree with degree distributions tt"','k^. If it"' and tt^ have finite first moment, then 



^ = \ ^ (1 _ max F"(t)), (8) 

' (/.«'(1) + /'(1)^ tern ^ ^ ' 

where F"- is defined by: 

Note that if (j)°'{x) = 4>^{x), we find the same limit as in Theorem [2J Note that it is not 
obvious from formula ([8]) that our expression for 7 is symmetric in a and b as it should. In 
the forthcoming Section [5.21 Equation (|30p gives an alternative symmetric formula for 7 which 
simplifies to ([8]) thanks to ([7]). 

Note also that our Theorem [3] computes the independence number of random bipartite 
graphs. Recall that a set of vertices in a graph G is said to be independent if no two of them 
are adjacent. The cardinality of any largest independent set of points in G is known as the 
independence number of G or the stability number of G and is denoted by a{G). By Konig's 
theorem, we know that for any bipartite graph G with vertex set F, a{G) + v{G) = \V\. The 
fact that a limit for exists, has been proved recently in for Erdos-Renyi and random 

regular graphs. The actual value for this limit is unknown except for Erdos-Renyi graphs with 
mean degree c < e. In this case, the leaf-removal algorithm allows to compute explicitly the 
limit which agrees with ([8]) with (p"'{x) = 4>^{x) = exp(cx — x). 

Motivated by some applications for Cuckoo Hashing [18], [15], recent results have been 
obtained in the particular case where 7r"(A;) = 1 for some A; > 3 and tt^ is a Poisson distribution 
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with parameter ak. These degree distributions arise if one consider a sequence of bipartite graphs 
with [amj nodes of type a (cahed the items) , m nodes of type h (cahed the locations) and each 
node of type a is connected with k nodes of type b chosen uniformly at random (corresponding 
to the assigned locations the item can be stored in). The result in this domain, obtained in |19j 
follows (see Section [5. 3p from our Theorem [31 namely: 

Corollary 4 Under the assumption of Theorem\^ and with 7r"(A;) = 1 for some k > 3 and vr^ 
is a Poisson distribution with parameter ak. Let ^ be the unique solution of the equation: 



and Oc = unSt\k~i ■ 



fc(l-e-«)'' 

for a < ac, all (except Op{n)) vertices of type a are covered, i.e. ^\8a \ > 1- 

for a > Uc, we have: 

#r — ^ 1 - - f^"'* + + f (1 - - ' (9) 

I I n^oo a \ k J 

where ^* = kax* and x* is the largest solution of x = (l — ^ . 



In words, Oc is the load threshold: if a < Oc, there is an assignment of the [amj items 
to a table with m locations that respects the choices of all items, whereas for a > Oc, such 
an assignement does not exist and ^ gives the maximal number of items assigned without 
collision. Note that results in [18], [15] are slightly different in the sense that for the specific 
sequence of random graphs described above (i.e. uniform hypergraphs), they show that for 
a < Oc all vertices of type a are covered with high probability. It is shown in [23] how to get 
such results from Corollary H] under the additional assumption that the sequence of graphs are 
uniform hypergraphs. 



3 The Monomer-Dimer model 

We start with the case of a finite graph G = {V,E). Consider a natural family of probability 
distributions on the set of matchings M(G), parameterized by a single parameter z > called the 
temperature (note that the standard temperature T in physics would correspond to z = e~^l'^ 
but this will not be important here): for any M G M(G), 

|y|-2|M| 

where Pq is the matching polynomial, Pciz) = Z]AfeM(G) statistical physics, this 

is called the monomer- dimer model at temperature z on G (see [21] for a complete treatment). 
We let Mq denote a random element of M(G) with law Hq, and we call it a Boltzmann random 
matching at temperature z on G. Note that the lowest degree coefficient of Pq is precisely 
the number of largest matchings on G. Therefore, Mq converges in law to a uniform largest 
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matching as the temperature z tends to zero. We define the root-exposure probability (REP) of 
the rooted graph {G, o) as 

TZz[Gi °] = /"g (° is exposed) . (11) 



3.1 Local recursions to compute 7 

Before starting with the proof, we explain (whithout proofs) how to compute 7 in ([3]). For 
a finite graph, our computation of 7 follows exactly the approach of Godsil [2D]. We recall 
Godsil's notion of the path-tree associated with a rooted graph G: if G is any rooted graph 
with root o, we define its path-tree Tq as the rooted tree whose vertex-set consists of all finite 
simple paths starting at the root o; whose edges are the pairs {P, P'} of the form P = vi . . .Vn, 
P' = vi . . . VnVn+i{n > 1); whose root is the single-vertex path o. By a finite simple path, we 
mean here a finite sequence of distinct vertices vi . . .Vn {n > 1) such that fifi+i G E for all 
1 < i < n. 

It is well-known since Godsil's result [20j that path-trees capture considerable information 
about matchings in general graph and are easier to work with than the graph itself. For a 
rooted graph [G, o], let T^q ^] be the associated path-tree and consider the corresponding system 
of equations (where u y v if u is a, child of v): 

Wv e T[G,o]' = 7^; (12) 

For any finite rooted graph [G, o], (jl2p has a unique solution in [0, IJ^I'^-"! and we denote the 
value taken at the root by Xo{G). Then Xo{G) is exactly the probability for the root o of being 
exposed in a uniform maximal matching. In particular, we have 

x.{G) 



KG) = 1 



This argument follows from [20] and will be a special case of our analysis. 

For an infinite graph with bounded degree, it turns out that it is not always possible to 



make sense of the local recursions (jl2p . However, our analysis will show that for any z > 0, the 
infinite set of equations: 

has a unique solution in [0, l]^l<3'°l and the value taken by the root is exactly TZz[G,o\ (which is 
the probability for the root o of being exposed in a Boltzmann random matching at temperature 
z when the graph G is finite) . Then our Theorem [11] will imply that for any sequence of finite 
graphs {Gn = {Vn-,En))n&^ satisfying \En\ = 0(114.1) and having p as a random weak limit, 

v{Gn) 1-Ep [lim,^o7^,] 

^ „ 1 



and \\mz-^QTZz[G, o] is actually the largest solution to (|12p . From a practical point of view, it is 
possible to compute an approximation of TZz[G, o] by looking at a sufficient large ball centered at 
the root o. Moreover our analysis will show that the quantity Ep \TZz\ is a good approximation 
of Ep [lim^_!.o 7^2] as soon as \En\ = 0{\Vn\) (see Lemma [T2]) . 



7 



3.2 Extension of the model on infinite graphs with bounded degree 



Let G — o be the graph obtained from G by removing its root o. Since the matchings of G that 
expose o are exactly the matchings of G — o, we have the identity 

^.[G,o] = ^^fi^, (13) 
Pg{z) 

which already shows that the REP is an analytic function of the temperature. The remarkable 
fact that its domain of analyticity contains the right complex half-plane 

IH+ = {z G C;5R(z) > 0} 

is a consequence of the powerful Heilmann-Lieb theorem |21^ Theorem 4.2] (see [14J for gener- 
alizations). The key to the study of the REP is the following elementary but fundamental local 
recursion : ^ 

n,[G,o] = z'\z' + y]n,[G-o,v]\ . (m) 



Clearly, this recursion determines uniquely the functional IZz on the class of finite rooted graphs, 
and may thus be viewed as an inductive definition of the REP. Remarkably enough, this alterna- 
tive characterization allows for a continuous extension to infinite graphs with bounded degree, 
even though the above recursion never ends. We let % denote the space of analytic functions on 
EI+, equipped with its usual topology of uniform convergence on compact sets. Our fundamental 
lemma is as follows : 

Theorem 5 (The fundamental local lemma) 

1. For every fixed z G IH+, the local recursion JT^ determines a unique TZz : Gd zM.+ . 

2. For every fixed [G, o] G Q^, z i— )• TZz [G, o] is analytic. 

3. The resulting mapping [G, o] G i — )• 7^(.')[G, o] ^T-L is continuous. 

This local lemma has strong implications for the monomer-dimer model, which we now list. The 
first one is the existence of an infinite volume limit for the Gibbs-Boltzmann distribution. 



Theorem 6 (Monomer-dimer model on infinite graphs) Consider a graph G £ Gd o-i^d a 
temperature z > 0. For any finite matching M of G, the cylinder-event marginals defined by 

2\M\ 

fihiMCM) = J] 7^,[G - W, . . .,vk-i},vk], 

k=l 

are consistent and independent of the ordering vi, . . . V2\m\ of the vertices spanned by M. They 
thus determine a unique probability distribution fj,Q over the matchings of G. It coincides with 
the former definition in the case where G is finite, and extends it continuously in the following 
sense : for any o G F and any sequence {[Gn,On])neN S Gd^ converging to [G, o], 

[Gn,On,Mn] [G,o,M], 
n— )-oo 

in the local weak sense for random networks, where 7W„ has law fJ,Q^ and M has law fiQ. 
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Although it is not our concern here, we obtain as a by-product the strong convergence of the 
logarithm of the matching polynomial, also called free energy in the monomer-dimer model : 

Corollary 7 Let {Gn)ne'N be a sequence of finite graphs with bounded degree admitting a random 
weak limit [G, o]. The following convergence holds in the analytic sense on IHI_|_, 



where Kp[TZs[G,o]] denotes the expectation under the measure p of the variable lZs[G,o\. 

A similar result was established in [21] for the lattice case, and in [6] under a restrictive large 
girth assumption. 

3.3 Proof of Theorem [5] : the fundamental lemma 

The local recursion ()14p involves mappings of the form : 



where d S N. In the following lemma, we gather a few elementary properties of this transforma- 
tion, which are immediate to check but will be of constant use throughout the paper. 

Lemma 8 (Elementary properties) For any d € N and z G ]H_|_, 

1- 'Pz^ maps analytically zM.j^ x . . . x zM.j^ into zM.j^ 

-2- \4'z,d\ "^s uniformly bounded by |z|/K(z) on zM^ x . . . x zM.^. 

From part 1, it follows that the REP of a finite rooted graph belongs to Ti, when viewed as a 
function of the temperature z. Part 2 and Montel's theorem guarantee that the family of all 
those REPs is relatively compact in H. Note that relative compactness also plays a central role 
in [21]. Combined with the following uniqueness property at high temperature, it will quickly 
lead to the proof of Theorem [5l 

The local recursion (|14p also involves graph transformations of the form [G, o] i— [G — 
o,v], where v ^ o. Starting from a given [G,o] £ Q^, we let Succ*[G,o] C denote the 
(denumerable) set of all rooted graphs that can be obtained by successively applying finitely 
many such transformations. 

Lemma 9 (Uniqueness at high temperature) Let [G, o] £ and z G H-^ such that 'R{z) > 






■+ 
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Proof. Set a = 2\z\/^{z) and /3 = ^(z) ^. From ()14p and part 2 of Lemma [8] it is clear that 
the absolute difference A = \TZl — 1Z1\ must satisfy 

A[G,o]<a and A[G, o] < /3 ^ A[G - o, 

In turn, each A[G — o,ti] appearing in the second upper-bound may be further expanded into 
f^Yliw^v w^o ^['^ — o — v,w]. Iterating this procedure k times, one obtains A[G, o] < j3^d^a. 
Taking the infimum over all k yields A[G, o] = 0, since the assumption lk{z) > \fd means 
precisely fid <\. □ 

Proof of Theorem O For clarity, we divide the proof in three parts : we first define a specific 
solution which satisfies (I14jl . We will then prove its unicity and check its continuity. This will 
prove parts 1-3 of Theorem [5j 

Analytic existence. Fix [G, o] g Q^, and consider an arbitrary collection of IHI_|_ — )• zM^ 
analytic functions z i— )• Tl^[H,i], indexed by the elements [H,i] £ Succ*[G, o]. For every n > 1, 
define recursively 

K[H,i]=z' (^z' + Y,n^,-'[H-zj}j , (15) 

for all z £ and [H,i] G Succ*[G,o]. By Lemma (H each sequence (z i— )• 7^"[iJ, i])„gj^ is 
relatively compact in Ti. Consequently, their joint collection as [H,i] varies in the denumerable 
set Succ*[G,o] is relatively compact in the product space '}{^^^'^ [G,o]^ Passing to the limit in 
([15]), we see that any pre- limit TZz- Succ*[G,o] — )■ zIH+ must automatically satisfy (fH|) for each 
z G IHI_|_. By Lemma m this determines uniquely the value of TZz[G,o] for z with sufficiently 
large real part, and hence everywhere in EI+ by analyticity. To sum up, we have just proved the 
following : for every [G, o] £ Q^^ the limit 

7^^[G,o] := lim 7e!^[G,o] (16) 



exists in 7^, satisfies the recursion and does not depend upon the choice of the initial 

condition 7^^ : Succ*[G, o] — )• zH-i- (provided that the latter is analytic in z G ^+)- 

Pointwise uniqueness. Let us now show that any S: Succ*[G,o] — )■ zHI+ satisfying the 
recursion (I14p at a fixed value z = zq £ EI+ must coincide with the z = zq specialization of the 
analytic solution constructed above. For each [-ff, £ Succ*[G,o], the constant initial function 

[-f^)^] •= ^iH^i] is trivially analytic from to zH^., so the iteration ()15p must converge to 
the analytic solution IZz- Since 7^"^ = S for all n G N, we obtain T^^^ = 5, as desired. 

Continuity. Finally, assume that ([G„,o])^^-^ g Qd^ converges locally to [G, o], and let us 
show that 

^.[G„,o] ^^7^,[G,o]. (17) 

It is routine that, up to rooted isomorphisms, G, Gi, G2, . . . may be represented on a common 
vertex set, in such a way that for each fixed A: G N, [G^, o]^ = [G, o]^ for all n > n^. By construc- 
tion, any simple path vi . . .v^ starting from the root in G is now also a simple path starting from 
the root in each G„, n > n^, so the 7^— valued sequence {z 1— )• TZz[Gn — {vi, . . . , 'yfc-i}i Wfc])„>„^, 
is well defined, and relatively compact (Lemma [8]). Again, the denumerable collection of all 
sequences obtained by letting the simple path vi . . .Vk vary in [G, o] is relatively compact for the 
product topology, and any pre- limit must by construction satisfy (|14p. By pointwise uniqueness, 
the convergence (I17|) must hold. □ 
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3.4 Proof of Theorem [6] : convergence of the Boltzmann distribution 



Consider an infinite [G,o] G Qd, and let i[Gn,°])n>i be a sequence of finite rooted connected 
graphs converging locally to [G, o]. As above, represent G, Gi, G2, . . . on a common vertex set, 
in such a way that for each A; G N, [Gn,o]k = [G, o],^. for all n > n/.. Now fix an arbitrary finite 
matching M in G, and denote by vi, . . . ,V2\m\ the vertices spanned by M, in any order. By 
construction, M is also a matching of G„ for large enough n. But the matchings of Gn that 
contain M are in 1 — 1 correspondence with the matchings of Gn — {vi, . . . , V2\m\}j and hence 

But [Gn — {vi, . . . , Vk-i},Vk] converges locally to [G — {vi, . . . , Vk-i},Vk], so by continuity of TZz, 

2M 

{M<ZM) y z-^^ r\n,[G-{vu.. ■,Vk-i},Vk]. 

k=l 

Proof of Corollary Analytic convergence of the free energy follows from Theorem [6] and 

Lebesgue dominated convergence Theorem, since for any finite graph G = {V,E) we have 



The uniform domination 



— Wz) provided by Lemma [H □ 



4 The zero-temperature limit 

Motivated by the asymptotic study of maximum matchings, we now let the temperature z — t- 0. 
4.1 The case of bounded degree 

We first use the results from the previous section to prove a version of Theorem [1] for graphs 
with bounded degree. 

Theorem 10 (The zero temperature limit in graphs with bounded degree) For any [G 

Qd, the zero temperature limit 

7^JG,o] = lim i7^JG,o] 

2->0 

exists. Moreover, TZ^ : Qd — ^ [0, 1] is the largest solution to the recursion 

-IN -1 

o-v,w]] I , (18) 



with the conventions = 00, 00 ^ = 0. When G is finite, 7^*[G, o] is the probability that o is 
exposed in a uniform maximum matching. 
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Proof. Fix [G, o] g Q^. First, we claim that z i— )• TZz[G,o] is non-decreasing on M4.. Indeed, 
this is obvious if G is reduced to o, since in that case the REP is simply 1. It then inductively 
extends to any finite graph [G, o], because iterating twice ([UD gives 



-1^ 



O — V,W\ 



(19) 



For the infinite case, [G, o] is the local limit of the sequence of finite truncations ([G, o]„)^gpj, 
where, for n > 1, [G, o]„ denotes the finite rooted subgraph induced by the vertices lying at 
graph-distance at most n from o. So by continuity of the REP, 7^z[G, o] = limn^ooT^z[G,o]n 
must be non-decreasing in z as well. This guarantees the existence of the [0, 1]— valued limit 

7^JG,o] = lim i7^JG,o]. 

Moreover, taking the z — t- limit in (jl9p guarantees the recursive formula (jlSp . 

Finally, consider 5*: Succ*[G,o] — [0,1] satisfying the recursion (fTSj) . Let us show by 
induction over n G N that for every [H,i] £ Succ*[G,o] and z > 0, 

S,[H,i]<n,[H,i]2n. (20) 
The statement is trivial when n = {TZz[H,i]Q = 1), and is preserved from n to n + 1 because 

T^z [H, i]2n+2 



\2n 



^ [^ + Y.[Y.S*[H-i-J,k] 



Letting n — )• c« and then z — )■ in (|20l) yields 5=,, < TZ* , which completes the proof 



□ 

This naturally raises the following question : may the zero temperature limit be interchanged 
with the infinite volume limit, as suggested by the diagram below ? 



n,[Gn, 



■n,[G, 



2-S-O 



■7^4G, 



Unfortunately, the recursion (jlSp may admit several distinct solutions, and this translates 
as follows : in the limit of zero temperature, correlation decay breaks for the monomer-dimer 
model, in the precise sense that the functional 7^* : Gd ~^ [0, 1] is no longer continuous with 
respect to local convergence. For example, one can easily construct an infinite rooted tree [T, o] 
with bounded degree such that 

lim i7^*[T, o]2„ 7^ lim t ^* [^i °]2n+i- 
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Indeed, consider the case of T being the graph on N rooted at = o, where two integers share 
an edge if they differ by 1. Then, a straightforward computation gives 7^*[T, o]2„ = 1/2 while 
T^*[T,o]2n+i = 0. Despite this lack of correlation decay, the interchange of limits turns out to 
be valid "on average", i.e. when looking at a uniformly chosen vertex o. 



Theorem 11 (The limiting matching number of bounded-degree graph sequences) Let 

p be a probability distribution over Qii. For any sequence of finite graphs {Gn = (Kd -£'n.))n6N 
satisfying \En\ = 0{\Vn\) and having p as a random weak limit, 

i^jGn) ^ 1 - Ep [7^,] 

\Vn\ n^oo 2 



In order to get our Theorem [H we need to remove the bounded degree assumption. This 
is done below. In the case where the limit p is a (two-level) Galton- Watson tree, the recursion 
(jlSp simplifies into a recursive distributional equation (RDE). The computations for these cases 
are done in Section O 

Proof of Theorem 

Let G = {V^ E) be a finite graph and M be any maximal matching of G. Then 
^ l{v is exposed in M) = \V\-2^ l(e G M). 

In particular, if p = V({G), we have the elementary identity 

IEp[^*] = l-^- (21) 
The proof of Theorem [TT] will easily follow from the following uniform control: 



Lemma 12 (Uniform continuity around the zero-temperature point) Let G = {V, E) 

be a finite graph. For any < z < 1, 



[n^] + <^f>[n,] [n,] . (22) 

\V\ log z 



Indeed, let p be a probability distribution on Qci, and let (G„ = (14, £'n))neN be a sequence 
of finite graphs with \Vn\ = 0{\En\), whose random weak limit is p. For each n G N, set 
Pn = hl{Gn)- With these notations, proving Theorem [TT] amounts to establish : 

Ep„ > Ep [7^*] . (23) 

71— >00 

However, since pn =^ p, and since each TZz, z > is continuous and bounded, we have for every 
z>0, 

n^oo 

Thus, setting C = sup„g|^ |^ and letting n — )• oo in (j22]) . we see that for any z < 1, 
Ep [TZ,] + < liminf E„„ [7^,] < limsupEp„ < E„ [TZ,] . 

log Z n^oo 
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Letting finally z — )• 0, we obtain exactly (j23p . and it only remains to show Lemma [T2l 
Proof of LemmaUM. Fix < z < 1. Since 2 i— )• Ep [TZz] is non-decreasing, we have 



Ep [7^*] < Ep [TZ,] < 



log z 



s-^Ep [TZs] ds. 



Use Ep [T^s] = |yj|;^(g^ to rewrite this as 



Ep [7^,] < Ep [7^,] < _ log -^"^^^^ 



|y|iogz "Pg(i)' 

Now, Pg(1) is the total number of matchings and is thus clearly at most 21-^1, while Pg{z) is at 
least ^l^l~2i^(G)_ Ugjj^g (j2ip . these two bounds yield to 



Ep [n.] < Ep [n,] < 



1 



V\ log 2; 



|y|Ep [71*] log z- log 2). 



This gives (122 



□ 



4.2 The case of unbounded degree 



In this section, we establish Theorem [T] in full generality, removing the restriction of bounded 
degree from Theorem II li To this end, we introduce the d— truncation {d G N) of a graph 
G = {y,E), obtained from G by isolating all vertices with degree more than d, i.e. removing 
any edge incident to them. This transformation is clearly continuous with respect to local 
convergence. Moreover, its effect on the matching number can be easily controlled : 



u{G'^) < v{G) < u{G'^) + #{f G V-AagGiv) > d}. 



(24) 



Now, consider a sequence of finite graphs (G„)„gN admitting a random weak limit {G,o). First, 
fixing d E N, we may apply Theorem [TT] to the sequence {G'^)nen to obtain : 



1-Ep, [7^, 



\Vn\ n^oo 2 



where pd is the d-truncation of p. Second, we may rewrite 



as 



i^iGi) v{Gn) 



IK 



Wr,. 



< 



#{v G Vn;dQgG^{v) > d] 



Letting n — )• 00, we obtain 



lim sup 

n— >oo 



1 - pd [n,] u{Gn) 



IK 



< p (deg(o) > d) , 



This last line is, by an elementary application of Cauchy criterion, enough to guarantee the 
convergence promised by Theorem [H i.e. 



y{Gn] 

\Vn\ 



-> 7, 



where 



r 1 - Ep, 
7 := lim — 

d— >oo 2 



(25) 
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Note that because of the possible absence of correlation decay, the largest solution TZ^, [G, o] 
is not a continuous function of (G, o) E Q. In particular, we do not know whether it is always 
the case that 

7 = i^, (2.) 

as established in Theorem [TT] for graphs with bounded degree. However, (|26p holds in the 
particular cases where we have an explicit formula for Ep[7^^,] which depends continuously upon 
the degree distribution as will be the case in Section [5l 



5 Computations on (hierarchal) Galton- Watson trees 
5.1 The case of Galton- Watson trees 

We now investigate the special case where the limiting random graph is a UGW tree T. Specif- 
ically, we fix a distribution vr G V{H) with finite support (we will relax this assumption in 
the sequel) and we consider a UGW tree T with degree distribution vr as defined in Section [2l 
The random matchings A^f^, z > are perfectly well-defined, and all the previously established 
results for graphs with bounded degree hold almost surely. However, the self-similar recursive 
structure of T gives to the fixed-point characterizations (fT4|) and (fT8]l a very special form that 
is worth making explicit. 

Before we start, let us insist on the fact that lZz[r\ (z > 0) is random : it is the quenched 
probability that the root is exposed at temperature z, given the random tree T. In light of 
Theorem [H it becomes important to ask for its distribution. Let V ([0, 1]) denote the space of 
Borel probability measures on [0, 1]. Given z > 0, £ V (N) and n G P ([0, 1]), we denote by 
&u,zitJ') the law of the [0, 1]— valued r.v. 

^ 

where M ^ v and Xi, X2, ... ~ /i, all of them being independent. This defines an operator Qy^z 
on V ([0, 1]). The corresponding fixed point equation = Qu,z{l^) belongs to the general class 
of recursive distributional equations, or RDE. Equivalently, it can be rewritten as 



where Xi,X2, ■ ■ ■ are i.i.d. copies of the unknown random variable X. Note that the same RDE 
appears in the analysis of the spectrum and rank of adjacency matrices of random graphs |llj . 
[12j . With this notations in hands, the infinite system of equations (jl4p defining TZz[T] clearly 
leads to the following distributional characterization: 

Lemma 13 For any z > 0, TZz[T] has distribution Qn,z{fJ'z), where is solution to the RDE 

IJ-z = Qn,z{l^z)- 

The same program can be carried out in the zero temperature limit. Specifically, given 
z^, v' € 7^(N) and fi £ V ([0, 1]), we define QuyifJ-) as the law of the [0, 1]— valued r.v. 

Y = f , (27) 
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where J\f ^ v, Mi ~ u\ and Xij ~ /x, all of them being independent. This defines an operator 
on V ([0, 1]) whose fixed points will play a crucial role in our study. Then, Theorem [TOl 
implies: 

Lemma 14 The random variable 7^*[T] has law 07r,5f(/^*), where /i^, is the largest solution to 
the RDE fi^: = 65f^5f(/i*). 

Recall that the mean of @n,Tf{f^*) gives precisely the asymptotic size of a maximum matching 
for any sequence of finite random graphs whose random weak limit is T (Theorem [TT]) . We will 
solve this RDE in the next section in the more general set-up of UHGW trees. Combined with 
Theorem [T] and a simple continuity argument to remove the bounded degree assumption, this 
will prove Theorem [2l 

5.2 The case of hierarchal Galton- Watson trees 

As in previous section, we first assume that both vr" and tt^ have a finite support. We can define 
a RDE but with some care about the types a and b. The corresponding results read as follows: 

Lemma 15 For any z > 0, conditionally on the root being of type b (resp. a), TZz[T] has 
distribution ©^62(^2) (^esp. Q^b where /i" is solution to the RDE: 

and fil = e^b^^ifi^^). 

For z = 0; conditionally on the root being of type a (resp. b), the random variable 7^*[T] has 
law 0^a 5fi,(/u") (resp. Q^b^-^a{^^^)), where /i" is the largest solution to the RDE 

^il = Q^a^M), (28) 
and /i^ is the largest solution to the RDE fj^ = Q^b 5fa(^^). 

We now analyze the RDE ([28]). We define: 

F^ix) = <^«(l-^''(l-x))-|^(l-/(l-x)-x/'(l-x).) (29) 

Observe that 

F^'ix) = ^/"(l - ^) (^'^(1 - - x)) - x) . 

Hence any x where F"" admits a local maximum must satisfy x = (/>'*(1 — (1 — x)). We define 
the historical records of F'^ as the set of x G [0, 1] such that x = (f>'^ {1 — 4^ {1 — x)) and for any 
< y < X, F"'{x) > F°'{y) (the latter condition being empty if x = 0). 

Theorem 16 If pi < ... < Pr are the locations of the historical records of F"", then the RDE 
[2^) admits exactly r solutions ; moreover, these solutions can be stochastically ordered, say 
. . . <^r, and for any i £ {1, . . . , r}. 



16 



• Mi((0, 1]) =Pi ; 

• Q^a^b{fJ,i) has mean F"'{pi). 

The proof of Theorem [16] relies on two lemmas. 

Lemma 17 The operators Q^a and Q^a are continuous (with respect to weak convergence) 
and strictly increasing (with respect to stochastic ordering) on ^([0,1]). 



Proof of Lemma 17. It follows directly from the fact that, for any n > and any rei, . . . , re„ > 
0, the mapping 



-1 



□ 



is continuous and increasing from [0, to [0, 1]. 

Lemma 18 For any fi £ V ([0, 1]), letting p = fj, ((0, 1]), we have 

2. if Q^a i^bip,) < fj., then the mean of Q^a ^b{fi) is at least F^{p). 

3. if Qi^a i^b{fj,) > fj,, then the mean of Q^a ^b{^) is at most F"'{p); 

In particular, if ^ is a fixed point of Q^a ^b, then p = (f)"- {1 — (j)^ {1 — p)) and Q^a ^b{^) has mean 
F^ip). 

Proof of Lemma[lR In equation ([27]) it is clear that F > if and only if for any i G {1, . . . , A/"}, 
there exists j G {1, . . . such that X^j > 0. With the notation introduced above, this 

rewrites: 

e^.,^.(M) ((0, 1]) = r(i-^\i-f, ((0, 1]))) , 

hence the first result follows. 

Now let X ~ /i, y ~ G5fa_5f6(^), A/"" ~ tt"-, N"- ~ vf'', and let S,Si,... have the distribution of 
the sum of a tt''— distributed number of i.i.d. copies of X, all these variables being independent. 
Observe that 

1^1=1 '-'i 



1 

c> -1 



1 











l...^f'',s^>o} 


-E 







J -'-{Vi=l...Ar<^,S,>0} 



'-{^i=i...Ar'^,Si>o} 



Then, @^a^^b{fi) has mean 
1 



E 



-AT" c -1 



1 + Ef=l S: 



k=l 



P(Vi = l...Af^,Si >0)-^A:7r^E 

cp-{l-^\l-p))-r'{l)E 
<A'^(l-^(l-p))-r(l)E 



S 



-1 



5-i + l + nJi 5, 



k-i „ _i l{5>0,Vi=l...fc-l,Si>0} 



5"^ + i + Ef=i^. 

Y 



.Sfa I {S>0,\fi=l...Sf'',Si>0} 



Y + S 



45>0} 
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where the second and last hnes follow from ([5]) and Y ~ Q^a 5f6(/i), respectively. Now, for any 
s > 0, X I— )• x/x + sis increasing and hence, depending on whether B^f^^f (;u) > /.i or Q^a ^b{fi) < fi, 
@^ai^b{fi) has mean at most/least: 



-0^(1 -p)) -0"'(1)E 



X 



x + s 



4S>o} 



with Xi are i.i.d. copies of X independent of ~ vf^ and M* = YliLi l{Xi>o}- Now if X' is 
the law of X conditioned on {X > 0}, and X- are i.i.d. copies of X' , by exchangeability, we find 



X 



E 



X 



x + T^^x^ 



^r— i{^*>i} 



pE 



X' 



X' + E -Li x[ 



Hence finally, depending on whether Q^if^[^) > fj, or Q^a i^b{fj,) < n, Q^a^i^b^fi), Q^a^^b{fi) has 
mean at most /least: 



b''{l-(p\l-p))-p^'"{l)E 



1 



I-{A^*>1} 



But using the definition ([5]) and the combinatorial identity (n + 1)(^) = (d+ one easily 

derive : 



1 



n>l d=l 



n\ p'^{l — p) 



n—d 



d d+l 



F'^{p). 



□ 



Proof of Theorem [76l Let p G [0,1] such that — 0^(1 — p)) = p, and define ^uq = 
Bernoulli(p). From Lemma [T8l we know that B^^a (/io) ((0, 1]) = p, and since Bernoulli(p) is 
the largest element of 7^([0, 1]) putting mass p on (0, 1], we have Qi^a i^b (hq) < hq. Immediately, 
Lemma [T71 guarantees that the limit 

fioo = lim \ e|„ (/^o) 

exists in 7'([0, 1]) and is a fixed point of ©^fa^fb. Moreover, by Fatou's lemma, the number 
Poo = Aioo ((0, 1]) must satisfy poo < P- But then the mean of Q^a i^b{fj.oo) must be both 



equal to F"-{p^) by Lemma [T8l with ^ = 



• at least F"'{p) since this holds for all Q^a^b o -.i,(/io), > 1 (Lemma [18] with /i 



We have just shown both F°'{p) < -F"(poo) and poo < p. From this, we now deduce the one-to- 
one correspondence between historical records of F^ and fixed points of Q:^a i^b. We treat each 
inclusion separately: 
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1. If admits an historical record at p, then clearly = p, so fioo is a fixed point satisfying 

/Xoo((0,l]) =P. 

2. Conversely, considering a fixed point n with /x((0, 1]) = p, we want to deduce that -F" 
admits an historical record at p. We first claim that // is the above defined limit Hoo- 
Indeed, ^ < Bernoulli(p) implies /i < /ioo (©s^a.^fft is increasing), and in particular p < poo- 
Therefore, p = poo and F"'{p) = F"(poo)- In other words, the two ordered distributions 
@^ai^b{fi) < 0^a5fb(/Xoo) share the same mean, hence are equal. This ensures fj, = fioo- 
Now, if q < p is any historical record location, we know from part 1 that 

i/oo = lim \ 0|a (Bernouni(q)) 

fc— >oo ' 

is a fixed point of Q^a^t satisfying i/oo((0, 1]) = q- But q < p, so Bernoulli(q) < 
Bernoulli(p), hence i/qo < /"oo- Moreover, this limit inequality is strict because z/qo ((0, 1]) = 
q < p = /ioo((0, 1]). Consequently, Q^a ^b{h'oo) < ^■k'^ n''i^^oo) and taking expectations, 
F"'{q) < F"'{p). Thus, admits an historical record at p. 



□ 

We may now finish the proof of Theorem [3l 

Proof of Theorem ; case of hounded degrees. We assume that Tr^ and vrf, have bounded 
support. Recall ([7|), so that we have A = ^sr^jpf^T^jy ; where A is the probability that the root 
is of type a. Theorems [11] and [16] and Lemma [15] give: 

v{Gn) ^ A(l - max,g[o,i] F"(x)) + (1 - A)(l - max^g^p^i] F\x)) 

Vj-i n— 5>oo J, 

where F°' is defined in (I29jl and is defined similarly by 

F\x) = /(l-0'^(l-x))-^(l-c/>"(l-x)-xr(l-x)). (31) 

For any x which is an historical record of F"", we define y = 0^(1 — x) so that 0"(1 — y) = x. 
Then we have: 

Hi-F^i^)) = x(i- ni - y) + r (1) { - ^'^^ -f}' - - vni - y)\ \ 




1 



{l-X){l-F\y)). 



By symmetry, this directly implies that A(l — max^.g[o^i] F°(x)) = (1 — A)(l — max^.g[o^i] F'^{x)) 
so that (I30p is equivalent to (j8]). This proves Theorems [2] and [3] for distributions with bounded 
support. □ 

Proof of Theorem [3 ; general case. To keep notation simple, we only prove Theorem [2] 
The following proof clearly extends to the case of UHGW trees. Let Gi, G2, . . . be finite random 
graphs whose local weak limit is a Gallon- Watson tree T, and assume that the degree distribution 
TT of T (with generating function (f)) has a finite mean : </>'(l) = '^Zn'^'^ri < 00. For any rooted 
graph G and any fixed integer d > 1, recall that is the graph obtained from G by deleting 
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all edges adjacent to a vertex v whenever deg(t') > d. Hence T'^ is a Galton- Watson tree whose 
degree distribution vr"^ is defined by 

Vi > 0, Vrf = TTili^d + h=0 ^ TTfc. 

k>d+l 

By Theorem (TJ Equation (I25p and our weaker version of Theorem [2] for distributions with 
bounded support, 

> lim min /(x), (32) 

ra->oo d-^ooxe[0,l] 

with </)rf(2;) = Ylk=o'^kx'' and 

/(x) = 1 - ^(1 - ^)'^'.(x) - \Mx) - l<Pd (l - ^ 

Also, as d — )• CO, we have (pd ^ (p and 0^ — )• 0' uniformly on [0, 1], so 

min g'^(x) > min g(x), (33) 

x€[0,l] n->oo ze[0,l] 



with g{x) = 1 — ^(1 — x)(l)'{x) — ^(j){x) — ^(p (^l — ^T^fy^ • Finally, combining and 
easily obtain the desired 

— j > mm (7(3; j. 

a;e[0,l] 



we 



□ 



5.3 Proof of Corollary |4] 

Note that in Corollary HI we divide v{Gn) by \Vn \ = \am\ instead of + = [amj + m, 

so that by Theorem [3l we have nwr > ™iiitefo il 1 ~ F"-{t). We have (/)"(x) = x^ , (j)'^{x) = 

^ak{x-i) gQ have: 



(0 



Let x* be defined as in Corollary [5] as the largest solution to 2; = (l — e '^"^)'^ ^. It is easy to 
check (see Section 6 in [23] for a more general analysis) that 

min 1 - F'^it) = min|l, 1 - F°-{x*)}. 
te[o,i] 

Setting ^* = kax*, we have |^ = (1 - e"^*)^"-*^, so that 

min l-F^'it) = mini 1,1 - - ( e"^* + Te"^* + —(1 - e"?*) - 1 | 1 . 
ie[o,i] [ aV ^ yj 

Since z 1— )• is increasing in z, we see that ^* > ^f and only if a > Oc and we get 



min 1 - F'^it) = 1 - l(a > aj- e"* + ^e"^ + -^(1 - e"^ ) - 1 • 
te[o,i] a V A; ' 
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Appendix : uniqueness and non-uniqueness at zero temperature 

Consider a sequence (G„ = {Vn, En))n£N of finite graphs whose local weak limit under uniform 
rooting is a UGW tree T. Let ^(t) = Yl T^nt^ be the generating function of the degree distribution 
TT of T, and for t G [0, 1] set 

Fif) = t<t>'{l -t) + <t>{l-t) + <t>[l- - 1- 
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From Theorem [T6l we know that there is a.s. a unique solution to the local recursion at tem- 
perature z = on T (correlation decay) if and only if the first local extremum of F is a global 
maximum. In that case, the convergence 

HGn) ^ 1 - max,g[o^i] F{t) ^^^^ 

I Vn I n-^oo 2 

can be obtained by a fairly standard compactness-uniqueness argument, without any need for a 
detour through the positive temperature regime. As it is not hard to check, a sufficient condition 
for the first local extremum of F to be a global maximum is that 0" is log-concave. This is in 
particular true in the Erdos-Renyi case, where (j){t) = exp(ct — 1), (c > 0). The corresponding 
function F is given in Figured] for various values of c. 




Figure 1: From left to right : plot of F for c = 2, c = e and c = 3. 

However, there are simple examples of degree distributions vr for which the function F has 
more than one historical record, implying the coexistence of multiple solutions to the local 
recursion (lack of correlation decay). In that case, a detour through the positive temperature 
regime is needed in order to establish the convergence ()34p . Here is a couple of examples. 




Figure 2: Plots of F for ^{t) = f + It^s (left) and (j){t) = §^t^ + ^t^o + ^t^oo (right). 
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